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ABSTRACT 


The nonlinear flap-lag coupled oscillation of torsionally 
rigid rotor blades in forward flight is examined using a set 
of consistently derived equations by the asymptotic expansion 
procedure of multiple time scales. The regions of stability and 
limit cycle oscillation are presented. The roles of parametric 
excitation, nonlinear oscillation, and forced excitation played 
in the response of the blade are determined. 
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I. Introduction 

The problem of the coupling between nonlinear oscillation, 
parametric excitation, and forced oscillation is of great interest 
in rotor dynamics . The present paper will restrict con- 
sideration to the flap-lag motion. Few authors have considered 
the nonlinear oscillation of a torsionally rigid blade. Young 
(Ref. 1) has used rough approximations to derive the solution 
of a set of nonlinear equations. Hohenemser and Heaton (Ref. 2) derived 
the nonlinear equations of motion in forward flight approximately 
and solved the equations by numerical integration. The present 
author (Ref. 3) used the equations derived by Hohenemser and applied 
the asymptotic expansion technique to establish the limit cycle 
oscillation amplitude and clarified the role played by nonlinear 
oscillation, parametric excitation, and forced oscillation of the 
rotor blades. In the present paper, the nonlinear 
motion of flap-lag in forward flight are carefully rederived, 
which is extended from the linear flap— lag equations in hover 
derived by Ormiston and Hodge (Ref. 4). A technique similar to 
that of Ref. 3 is applied to treat these new equations. Extensive 
results on the forced response, stability boundary > and limit 
cycle amplitude will be presented. 

II. Formulation 

In deriving the equation of motion, we shall make the follow- 
ing assumptions! (1) the pitch angle, 0, the static equilibrium 
coning angle and the static equilibrium lag angle are small, 

so that 0 2 , 0 2 and £ 2 can be neglected as compared to one; (2) the 
effective angle of attack is small and the two-dimensional quasi- 

If — ' 

All of the symbols used are defined in Appendices A and B. 
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steady aerodynamic theory can be used; ( 3 ) the ratio of the 

profile drag coefficient to the lift slope is small and can be 

neglected compared with unity; and (4) nonlinear terms up to the third 

order will be retained. However, since the aerodynamic damping is much 

higher in flap than that in lag, the perturbed motion in lag is much 

larger than that in flap; thus, it is expected that the nonlinear 

terms of the perturbed lag motion 5 will be more important than the 

perturbed flap motion 3. It is also expected that the nonlinear terms 

are more important in the lag equation than in the flap equation.* Thus, 

2 2 

the nonlinear terms such as 8 C, 8 C , 6 ? and their time derivatives 

2 2 

will be retained in the flap equation while 6 / 8 S, 08 S etc. will 
be retained in the lag equations. 

From the above assumptions, it has been shown (Appendix A) 
that the equations of motion in flap and in lag are* respectively: 

& + %, * + P - if 

4 s ^T e 5 + 3^- 

(1) 

r 1 ■■ — > 

This is indicated by the results of ASRL 166-1 (Ref. 3). 
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■+ ^ 5-5 •+ + P 1 X P 

= a.p£ ■*{»&,-$«>) P>t •+ -apii 

-►i 4 f- (4 a. •*ce)sM'-. (4P.&-2M w>^] + x u{[^.iC)« 

3 L ^ CC ^ 

_ i.^!>s-5e0>i.'['+[(ic-|e)|i*-|f ! 'o(2p-9i)- c es]‘ <>t /']j 

-* l (f&- JC )P 5 *f 

+ [j| (ip, - ^| ©C fit + &S>] (2) 

-v xy |* pi «>*f - (*.& ^ *+- 2 % 


• J 

where ( ) = ip = St, y is the Lock number, y is the advance 

ratio and the rest of the parameters are defined in Eqs . (A.9), 

(A. 22), (A. 23), and (A .24) . 

The critical condition of the system is determined by the 
existence of a small amplitude undamped oscillation for the linearized 
version of Eqs. (1) and (2). Namely, 

(*> P ^ “ *2*5 " Q 

f - 1 \ i * -»s * t jip (b - Y (i = o 


( 3 ) 



in (3) , we have the characteristic equation for p 


where i = ^ - 1 . The critical condition corresponds to p = f2 F which 
is called the flutter frequency. Substituting into (4) and using 
subscript "c" to denote the critical condition for the parameter, 
one obtains 

-4 -+35cCi+.fip )+J? (Ke fe c +x)]/ <5) 


{6) 

* 

If all the physical parameters such as 

and are given, one can determine the critical pitch angle, 

say 0_ from (5) and (6) . It can be shown that for 0<0 all the 
c c 

small amplitude motion will be damped and for 0>0 c there exists 

a branch of solution of (3) , which will grow exponentially. 

Ormiston and Hodge, Ref . 4, have shown that when e, and 

2 

^5$ =0, it is necessary to have 0<fi^<l for the existence of 0 c « 
Indeed in this case, it can be shown that 


i > -ftp 


(7) 


^ he e ^ ast ^- c coupling effect, no expression for it is derived 

at present. It will be treated in a forthcoming report. It was 
shown in Ref. 4 to have an important influence on flap-lag insta- 
bility. 
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We define the differential operators L and M, the functions 
f's and h's, and various quantities a's, 8 . . , £ • • , N, P's and <'s 

1 J 1 J /V 

in Appendix B. All these quantities are to be used in the follow- 
ing sections for the perturbation solutions. 

Ill, Perturbation Solution below the Critical Condition 

If the angle of the pitch is less than 0 , for a given set of 
physical parameters, all the solution of (3) are damped. A per- 
turbation solution of (1) and (2) can be established by assuming 

( 8 ) 

$ = -* k — 


Substituting (8) into (1) and (2) , requiring all terms of the same 


order in y satisfying (1) and (2) separately, we obtain 

0 (y) 

U^,) =■ -f.K 


L Cf’i ^k') - + 

M( p.,5,) - (% M + K 5 (t 8; p, -5, ) + A) 


(9) 


(10) 


where L, M, f's and h's are defined by (B.l) through (B.17) in 
Appendix B. 

Since the homogeneous solutions in (9) , (10) , etc. are damped 

(below critical condition) , we only have to consider the particular 
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solutions . Using (B.4) and (B.5), we may write (9) in the form 


t . T>* W 


L (6. , ~ -*• Vn t 


M Pge 1 '*' -v? s f e 


-ilk 


( 11 ) 


where P ? and P g are defined in (B.34) and (B.35). ( )* denotes 

the complex conjugate of the corresponding quantities. The par- 
ticular solution for B 1 and ^ is simply 

* ■ ' < ^ ( .-f 


($)- 


> %7. 


(12) 


where 8 ^ and ^\2 are de fi ned in (B.36). Substituting (12) into (10) 
and using (B.5) through (B.10), we have 

it 0+ 'iit 
-v -v r.s 

(13) 




and the particular solution for 3 and £ is 

* 4 



where p i 3 »*»* p 16 ^ ^25 f ? 25' 6 26' and ? 26 are defined in (B.41) 
through (B. 48) in which Z Q Q , $ and f are set to 1. 

By a similar procedure one can obtain the solutions for 
higher order expansions. 

IV. Pertu rbation Solution in the Neighborhood of the Cr itical 
Condition 

In the previous section, we only seek the particular solution 
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of 3 and If the pitch angle is larger than 0 C , at least one 

branch of the homogeneous solution of 

=o 

MtfcS) 

is undamped. This undamped solution must be included in the solu- 
tion of S-l etc. In the neighborhood of 0 C , the value of $ 12 ' 

5 8 , etc. becomes very large, i.e., the solution in (12) or 

M J 

(14) etc. are large (see B.36 and B.45 through B.48), the expan- 
sion in (8) will break down. Under these circumstances, we 
shall establish a different expansion by perturbing the system 
near the critical condition. Introduce a small parameter e and let 

for the dependent variables and let 

9 = + - ' 

= v + - 

ft a (16) 

\ = V 4 + •• * kc 4 •* • • • 

9 yf ■%«.* r + i f ' 
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The subscript "c" refers to the value at the critical condition. 

We also introduce different time scales. 

^ f ^ ^ ; w\ -=■ o v i , “ ■ (17) 

then 

■* % ■ ■* ■ 

- it, >r» 

X * 

■=- 4 AS jT 4 

For different values of Qp, different order of magnitude of 

parametric excitations and forcing functions, i.e., different order of 

2 2 . . 
y, yf^, yh^, y h 2 and y h 2 , are permitted in the expansion series 

of (15) . We write 

7 *# 

/ 1 2 ct 0>9 ) 4>o,e) 

Assuming f£, h^, f£ and h£ are functions at most of order one, 

we expand y, 9£ and ^ in the following forms: 

1 

yU r. U, 4 ' 

a£ = *>„ ■* 

1 - 


( 19 ) 
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The constraint in the choice of the value of the coefficients in 
(19) is to make the expansion series (15) a periodic solution 
in the fast time variable Tq. Since we shall only consider the 
solution up to order e ' , we may, in general, have only one of 

the coefficients as nonzero. In particular, within the above 
mentioned constraint, only the leading coefficient is chosen to 
be nonzero. 


Substituting (15) through (19) into (1) and (2) , requiring 
terms of the same order in e satisfying (1) and (2) separately, 
we will obtain the differential equations for the determination of (3^, 


5 etc. The equations can be expressed in the following form: 

0 < e 1/2 > 


0(e) 


M ( £>&) " />**) (21) 

- tMS.) 

' 2 ) 

+ + l\ 7 c^, ; \ • a ^ 4... 


where L, M, f's and h's are defined in (B.l) through (B.17). In 
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this section, the parameters used in L, M and N [(B.l), (B.2) and 

(B. 18) ] are referred to the value at the critical conditions. 

The homogeneous solution of (20) Is the same as (3). At the 
critical condition, there is only one branch of the solution 
which is undamped. Therefore, we may just include that particular 
branch in the solution of 6^ and and 



where 3^ and are defined in (B.30), (B.31), (B.36), (B.47) 

and (B.48). In (23), the y's, %'s and ' s are so chosen such 
that there are no secular terms in 3's and £'s, i.e., the solu- 
tions are periodic in the variable t q . A substitution of (23) 
into (22) yields 
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= ( e* +K p„) a A A* ] e t%r ° 

+/', [§£ P 7 + Pjs AA* '•/»? P 27 ] e lT » 

+ >!>.?<, A*e u, '- RF ' )r % XP» /yV (z ' J "^ r ° 




■*«&,« * e“ r °+ 

Mtp„\,)= [P^ - te»?»^R*)A -TtA'A^e 1 ^ 


(24) 


■f 




it* 


■+ M»P. A*e + A* A* e l fa 

■+ a>,? jo A^e lC '- iSl P r » + / u, f^A’e 1 fa5l f'^ T '’ 

+ + rt r ^P».« ,;c * + 

In the right hand side of (24) , many terms have the same 
frequencies as Q„, i.e., the solution of 3 3 , s 3 will have secular 
terms, i.e., 3 3 , £ 3 ~TQe lfi F T 0, unless the right hand side is 
orthogonal to the solution of the adjoint operator 


m 


-iBpt; o 


of L and M (Ref. 3) , i.e. 


\{* LC^') + H t a * v ' <tr 0 - o 


(25) 


in which v is defined in (B.63). Equation (25) will provide us 
the equation for A which can be functions of t 2 , i 4 ,...etc. 
Consider the following cases: 

(1) Hovering y(=y.=y 2 )= 0. In this case, all the forcing 


functions and the parametric excitations (terms with coefficient 
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y 1 or y 2 in (20) - (25)] are zero. To suppress the secular 
terms of an ^ (25) is equivalent to requiring 

Vr(?,^r T A ^5 f\(f) ' ?* T5c % + A+Pfc A"A V ^ O 


or 


3A a 
or. 


e-fci A ■+ 


(26) 


where v, < 2 , < 3 are defined in (B.63) through (B.65). Equation 
(26) is in the same form as that of Eq. (49) in Reference 3. It 
can be solved easily. The solution is 

i+ 


A = e& 


(27) 


where ( ) R and ( ) ^ are used to denote the real and jthe imaginary part, 


+-*. * ssM* 5/- 


(28) 


and p n and tf> n are the values of p and 4> at t o = 0. The conclusions 


concerning the stability of the solution is the same as that in 
Sec. IVa of Ref. 3; namely, (a) if 0 2 <2 R >O (above critical condil 
and k 3r > 0, the blade is unstable; (b) if 0 2 < 2R <O (below critical 
condition) and k 3r <0 the blade is stable; (c) if 0 o k ot ,>O and 


2 2R 


K 3R <0 ' there exists a limit cycle oscillation with amplitude being 




JL 

2 . 


( 29 ) 
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(d) if 9 o k ot ,<0 and k >0, there exists an unstable limit cycle. 

2 2R jK 

The blade is stable if the disturbance in | A | is smaller than 

that in (29), and the blade becomes unstable if the disturbance in |a 

is larger than that in (29). Let 0 2 =1, ky (15), (16) and (22), 

1/2 

the amplitude of the blade motion is (0-0 c ) A. 

It should be noted that k 2r is positive for the present 
problem. Therefore, the origin (A~0) is unstable above the 
critical condition (0 2 >O) , and is stable below the critical 
conditions (0 2 <O) . 

(2) y^O, ft F =l/2+ev, where v is of order one or smaller. 

In this case 3 23 and C 23 in the second equation of (23) can be 
very large or even singular, i.e., the solution of 0 2 and ? 2 will 
have secular terms. This is because the nonlinear coupling of 
the first order homogeneous solution, 8^ and with the forced 
response due to f.^ and h^ and its coupling with the parametric 
excitation are in resonance with the system. The difficulty can 
be avoided by requiring 

/t.-O (30a) 

and the other coefficients in (19) may be taken as 

i 

32, *se x = % - f o 


i.e. 


x^^\ 


(30b) 
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In (24), terms proportional to e^F T 0 and e"’" T 0 cause 


the secular terms in g 3 and x 3 . To suppress it, it is required 
that (25) be satisfied which is equivalent to requiring 

’dfS 


* W, + A h*e zi * Zl -PsAV] 




+ ^2. *t 2 + A* - 0 




or 


+ A ytjA“A* 


(31) 


Y)x 

changing variable by A=Ae 2, it becomes 




(e 2 v^u>)A -AiH 5 a*4 K.Va* 


(31) 


where k's are defined in (B.65). Equation (31) is in the same 
form as (10) of Ref. 3. Thus its behavior can be investigated 
in the same manner, that is, if 


A < I 


(32) 


(33) 


the origin, | A | =0 is a focal point; if 

I K zitWl , _ Ika^+cvl 

IM < ^ i"k 5 1 

the origin is a node; if 
| Ki +i^ 1 

— Hm ^ 

the origin is a saddle point. Since (31) is homogeneous, if 

* 'IfJlK.l*- ( n 21 e l j < o 

the solution near the origin (A=0) is stable. In this case, if k 3r <0, 


(34) 
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there will be no large amplitude flutter response, i.e. the blade 

will simply respond to the forcing function (no flutter mode) due 

to advancing flight. If e 2 <2R K 3R > 0 no close tra J ectory for A can 

be formed on the phase plane. In this case, if ^®2 K 2R <9 ^ ' 

the solution is expected to converge to one of the singular points 

of (31); if k 3 R >0 ( 0 2 k 2R >0 ^ ' the solution is unstable. If 

8 k k_ <0, a closed trajectory for A can exist on the phase 
2 2R 3R 

plane. The maximum distance of the trajectory from the origin 
is bounded by 




Ksr 



KinOt . 
*3R I 


(35) 


In this case, it is expected that if k 3r <0 ( 0 2 K 2R >0 ^ ' the closed 

trajectory is stable and if k 3r >0 (® 2 K 2 R <9 ^ ' ^ uns ^ a ^^ e * 

Let y 2 =l, i,e * e=y ' e 2 =(e- 0 c ) / y and T 2 =y1 ^ * By (15) ' (23) 

1/2 

and (31), the amplitude of the blade motion will be y A(0-6 c /y, y^) . 

For fixed ( 8-8 ) /y and as y increases, the amplitude will grow as 
c 

y 1 / 2 and will vary more rapidly with time scale ~l/y. A will tend 


9 ~ 6 c k 2R 
y *■>» 


as (e-8 c )/y^<», if e 2 K 2R >0 and K 3R <0 * 


(3) y^O, f2p=l+ev, where v is of order one or smaller. In 
this case, B ±2 and ? 12 in the first equation (23) can be very large 
or even singular, i.e., the solution of 6 ^ and 5 ^ has secular 
terms in t q . This is because the forcing function and h 1 have 
almost the same frequency as the flutter frequency. This diffi- 
culty can be avoided by requiring both f^ and h^ to be small; 
namely, ^=e or 
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ae 0 - ae, - o 
ae* 


(36a) 


in (19). This is a reasonable assumption because both and h^ 
are proportional to the pitch angle which are indeed relatively 
small. We may take the other coefficients in (19) as 


?.-l 


(36b) 


i.e. 

ae =■ * 

-I 

(In reference 3, y was taken to be y^e 2 ^ 2 , to be one. The 
present expansion is valid for a larger range of y.) 

In (24), the terms proportional to e-^^F T 0, e-^ T 0, e-'*' T 0 
and e-^ 2fl F 1 ^ T 0 will cause the secular solution for 3 3 and 
Equation (25) , in this case, is equivalent to requiring 

» [?,$, />>?«') A *t> £ ftV v * 

2* 

+/>. (T« kh tj i ^ v> p >7^ e " £vt * + /‘^iAV 2 ‘ 1,r *J 
j-y. ^ 12b) e' iori -*-/*? 1«A* e' 1 “ ITl = o 


or 


H = c w Z & 2 vX) A + K 1 ft + * ,(k lo A A*r ju>, ,4 *f)e 

-+ vi* a" a* -v /a 4 K,i A 1 " e L 




( 37 ) 



-17- 


Equation (37) is quite complicated. The complex amplitude 
"A" depends on y^ 0 2 , and v which are coefficients in the expansion 
of the parameters. Their physical meaning is not quite obvious 
because £ itself is just a conveniently chosen small parameter. 

To see the effects of pitch angle, advance ratio and the value 
of the flutter frequency, it is better to rewrite (37) in terms of 
physically more meaningful quantities by defining: 


A = Ae 




t z A 


e CCn,-OY 


(38) 


where 2 | A | is the amplitude of the first order solution of the 
lag motion in radians. Substituting (38) into (37), and using 
(16) - (19), we have 

|| = [ +c£VD] A + A* Ut A + M* 3 


The third bracket is the contribution due to the nonlinear 
coupling of the flutter mode and the parametric excitation and 
the term < 3 A 2 A* is from the nonlinear coupling of the flutter 
mode. If the nonlinear terms are neglected, we obtain the result 
for the linearized equations of (1) and (2) , i.e. 


A 

4^ 


2,. T * 


~t u Hu 

r . . _ 2 >* 
The homogeneous solution of (40) is in the form A Q e * 


( 40 ) 
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\a.= J V'*J l ‘Xxl 1 ‘ 

Therefore, the linearized solution, i.e., the origin (A=0) , is 
unstable if RefX^ or Re(X 2 )>0. Its stability boundary is given 
by 

Rfi. [ “ 0 (41) 

Since (39) is inhomogeneous, regardless of the range of the parameters, 
the solution for A will not be zero, i.e., the flutter mode is 
always excited. 

To examine the behavior of a nonlinear solution, one must use the 

full equation (39) or (37) to investigate the characteristic of its 

singular points. One can then estimate its behavior by the method 

used in Ref. 3. It can be shown that if y is sufficiently small 

and (9-6 )/y 2/ ^ 2 is of order one or larger, it reduces to (73) of 
c 

reference 3, i.e., |A|r/y ^/ 2 for small value of y and the time 
scale for the growth is proportional to l/y 2//2 . For a larger value 
of y, i.e., (0-0_)/y^ 2 is of order one, |a| will be of order 
y and the time scale for growth is proportional to 1 /y. 

(4) y^O, fip=2+ev. In this case, the forcing functions f 2 
and h 2 and the nonlinear coupling of the forced response due to 
f^ and h^, and its coupling with the parametric excitation will 
cause secular terms for $ 2 and ? 2 , i.e. 625 and ?25 can l ar 9 e 
or even singular in the second equation of (23) . The difficulty 
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is avoided by setting 

ae* =■ To = ° 
<*» 


(42a) 


The other parameters may be taken as 


at**'*,*® 

We have ^ 

/** /\£ l 

£^ 


(42b) 


(43) 


X T 

In the right hand side of (24), only terms proportional to e F 0 
and e 2 ^" T 0 will cause secular solution of 3 3 and £. 3 . To suppress 
the secular terms, (25) is equivalent to requiring 



or 


K 


IZ 


^ =(jMi+jA a *OA-+ *3^* +K-^e 


.ivt^ 


(44) 


Transforming back to the physically meaningful quantities by 


A -^A 


vt 2 


=. I 


1 


C%-2)(P 


we have 


( 45 ) 
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The behavior of A near the origin (A=0) depends only on the sign 

2 

of Re [< 2 ( 0— © c ) +y <g], i .e., stable if it is negative and unstable 
if it is positive. One can examine the general behavior of A by 
the method used in Sec. IV of Ref. 3. The results are the following 
Let 


If ^R < 3R >0f there exists no closed trajectory for A in the phase 
plane. Since the sign of 3 |a|/9t 2 is the same as < 3r for large 
values .of | A | . If < 3r <0 (<S 2 r <0) • the solution will converge to 
the singular points of (44). If < 3r >0 (6 2 r > 0) , the solution is 
unstable. If ^R^R^' there exists no trajectory for A on the 

•'N 

phase plane in the region | a| <- 6 2 r / 2< 3r . However, in this case, 
it can be shown that the maximum value of the closed trajectory 
from the origin is bounded by 


» C) i ML, * i 


(46) 


where Pq(> 0) and p^(<0) are the largest and the smallest real 
roots of 


1M N*o 


(47) 


It is expected that if i< 3r <0 (^X)) • the solution converges to 
the closed trajectory and vice versa. 

0 _ 0 

For fixed 0-0 and sufficiently small y, say -^- 7 ^ 0(1), or 

0-e c _ y 

~ >:> lf equation (45) can be solved similar to (26), with terms 

2 . _ 
proportional to y being neglected. In this case, A is of 
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order (Q-Q for y 1//2 , which is the same result as obtained in 
v> c 

Ref. 3. For larger value^of y, say (8-0 c ) /y~0 (1) the full 
equation (45), must be considered. The solution of A is of 

order y. 

Since (45) is inhomogeneous, A will not be zero, i.e.. 


flutter mode will always be excited. 

(5) y^O, ft F =l/2+ev,l+ev or 2+ev. In this case, there is 

no secular term for the solution of (8 1 , ? x ) and (6 2 ' ^ 2 ) * 
Therefore, we may take 


as. = - • 


(48) 


a * = h\ £Z 


(49) 


To suppress the secular terms of B 3 and c 3 , there are three cases 
(a) ft p =l/3+ev. In (24), the terms proportional to e ± lQ F T 0 
an( j el 1 (l _ 2Sp) T o _ will cause secular terms. Therefore, (25) is 


equivalent to requiring 

v[y, (e.T i +KPa')A*?6 a*V 3iVI * 3 


or 


ovz 


( 50 ) 
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Near the origin (A=0) , the solution is 


+/A i KfeR <• O 


unstable 

stable 


( 51 ) 


The general behavior can be examined by the method used in 
Ref. 3. Let 


S - K >2 “ 3 } 


(52) 


If 6 r k 3r > 0, there exists no closed trajectory for A in the phase 
plane. If S R < 3R <0, the maximum distance of the closed trajectory 
from the origin is bounded by 

[ I- S* I - Hi 1 . 

2 ) N ^ | KjR I lG,p| J x lyh^ 


< 


* A j ^3^ K5R 


j \ 


(53) 


(b) ft R =3+ev. In (24), the terms proportional to e-^F x 0 
3 i t 

and e 0 will cause secular terms. Therefore (25) is equivalent 
to requiring 


* A + p*. Ti = 0 

or 

^ - (H»©x -r M? A + A J ft* + 


( 54 ) 
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Equation (54) is similar to (45) , it can be examined accordingly. 

(c) Others: In (24), only the terms proportional to 

e 1 ^F T 0 will cause secular terms. Therefore, (25) is equivalent 
to requiring 

— (55) 

Equation (48) can be solved in a manner similar to that used 

for (26); therefore, the conclusions in Sec. IV. 1 hold for (55) with 

2 

< 2 0 being replaced by < 2 ®2 + y K 6‘ 

V. Uniformly Valid Expansion 

The solutions given in Secs. Ill and IV can be combined to 
establish a single expansion which is valid for small advance 
ratios p and all range of pitch angles 0 up to a small neighborhood 
above the critical conditions. We shall first show that there exists 
a range of 0 where the expansion given in both Secs. Ill and IV are 
valid, i.e., to establish the matching of the solutions. Then 
using a technique similar to the one used in singular perturbation 
theory (Refs. 5,6), a uniformly valid expansion can be constructed. 

Let us use subscripts "b" and "n" for the solution in 
Secs. Ill and IV to denote respectively the asymptotic expansion 
for 0 below and near the critical pitch angle 0 c , i.e. 



( 56 ) 
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the relation between e and y will depend on the cases (see Sec. IV) 
considered. Let e be a small parameter such that 


°< l , 


« 


(57) 


We shall show that for 0<0 and 

c 


e - e t = - f 

both expansions in (56) match, that is. 




o 


as £,e->-0 with t 2 being fixed and a=l/2, 
valid expansion is simply 


15 ) 


/ 


f 


P 




(58) 


(59) 


1,... Then the uniformly 


P 

S 7 C 


(60) 


The last term of (60) is the common part in the expan- 
sions (^) and (^) K . 

£ n C b 

We shall consider the four cases of Sec. IV separately. 
(1) Hovering, y=0. In this case 


•P 

5 



o 


S 

and fl) n is defined in Sec. IV. 1, and is proportional to A which 

is the solution of (26). Since k 5 0 o is linearly proportional to '0 ? 
0-0 . ^ 2 
(=—— £ =-•*■) (see B. 59 , B.29, B.39, 16 and A. 24), 

its real part is a large negative quantity. The solution of (26) 
can be approximately expressed in the form 



( 61 ) 


A = A.e Mfl 'A.e K,S *^ 


* 3 

where A Q is a constant. A is exponentially decaying, so is ( ? ) n 

i.e., (59) is automatically satisfied. We should point out that for 

the conditionally unstable situation of Sec.IV.l.d the approximate 
solution (61) may be invalid for a sufficiently large initial dis- 
turbance. However, this situation may be avoided by having sufficiently 
small £, thus seriously limiting the range of validity of the uniformly 
valid expansion. 

(2) n F =l+ev, (Sec. IV. 2). The expansion (®) ft is in the 
form similar (up to the order considered) to that of the previous 
case except that A is now governed by (37,) . In (37) , within the 
range given by (58) and y^ of order one or smaller, the real 

r 

part of K 2 0 2 (~-l/e) is a large negative quantity; the solution for 
A can be approximately expressed in the form 


A * A e 


-L 


(62) 


[This is valid only if 


2 1 2 
y |k 3 k 4 


i(n F -D+K 2 (e-e c )| 


!<<!•] 


where Aq is a constant to be determined from initial conditions. 
The expansion near the critical value 0 C can be written as 


/ - 0 \ >» K<|. /g\ /{L\ 

( l \ I Kz(e ®c)+ \ I ) **l Ki( 9 - 8 c) n 


X 




Ki 


Ki^-©c)+Ua r o 
4 0 (t*) | 






( 63 ) 


A q depends on the initial conditions. 
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We shall examine the expression (^), by expressing 

i = JTp - tV 

where v is a constant of order 1. Expanding the parameter in the 
form of (16), we have, (B.18, B.19 and B.36) 

3X0 = - l (-flp+i -4 -k 'P i ir»1^j+ oil) (64) 

^Vnu>^= 4 --J /PA 

*12 'b ^ / 


IHI 'i 


\r 


j / # \ 

\ /UJ 


-ISl^jLS [f^\. A/J_ _f/u\ ( ® ! 

Cn^Xp^^l' l ' 'j ~ i(K z 0 a 4-iv) 1 1 * h °^j 


From (14), (B. 47) , (B.48) 


rV T, 


( 66 ) 


= i(s f -'>K,(6-eol (^)„ + 

\ n i / ^ \ « — » J J 1 _ 


( 67 ) 


S2Z. # ” 

o 

Then the expansion (^) fa can be written as 


<*}<*? 
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( 1 ^ ~2 &. 

\3 


C(a r O^'&<) 


o(“K 


i + 


1 2 ^ 
M 



a* 

2 


_J54 3 / (V 

^ [6-8^1 ^ ^ ai j 




0(£‘) 


(68) 


Comparing (63) and ( 68 ) # evidently the matching requirement (59) 
is satisfied. The common parts of the two expansions within the 
order of approximation in the expansion are the same as in ( 68 ). 
Therefore, we have the uniformly valid expansion 



and 3 25 ^ C 2 s • ^26 and ? 26 are defined in the same wa Y as 3 25 

etc. in (B.47) and (B.48), except that and C ^2 i s replaced 

by r i2 and C 12 and -I . 

(3) n =2+ev. Similar to the previous subsection, it can 
F 

be shown that 8 and £ can be expressed in the same form as (69) 
in which 



where 
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(fc)' (%) 

with se 0 = l 5 0 =JE='f=i in (B. 41) - (B. 44 ) . 

For the other cases, the expansion near the critical con- 
ditions matches that below the critical condition, i.e., the 
uniformly valid expansion is of the same form as that in Sec. IV 
except that when evaluating $ 12 , c 12 , 6 25 , ? 25 , 3 26 and ?2g , 
the parameters in the matrix N (B.18) are referred to their value 
at the corresponding pitch angle 0, rather than using the value 
at 0 as done in Sec. IV. 

L* 

VI. Conclusions and Numerical Results 

To summarize the results for the pitch angle 0 in the neighbor- 
hood of 0 : 
c 

(1) The value of < 3r , defined in (B.65), characterizes the 

behavior of the blade at large amplitude motion because the sign 
B I A. I 

the same as that of when |a| is sufficiently 
large. Therefore, if k 3r has large positive value, the nonlinear 
effect is strongly destabilizing and vice versa. It is a very 
undesirable situation to have the blade with as positive, 

OX\ 

because instability can always be excited if the disturbance is 
sufficiently large. In practice, large disturbances are likely 
to occur due to gust or maneuver. The value of < 3R is independent 
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of the advance ratio; however, it depends on 0, ^ ^ t 

3 in a complicated way. In the two cases considered in 

p • c • 

Figs. (2) and (3), when 9^> 1 +0,*, < 3R is positive or has a very 

small negative value, i*e., nonlinear effect is destabilizing or 

is very weak. In the region < 3R generally negative, 

2 2 

and has a larger value than that for £2^ > 1 + i.e., the nonlinear 

effect is stabilizing and is stronger. 

(2) When flutter occurs, the amplitude of lag motion is 
much larger than that of flap. This is due to the fact that the 
total damping in lag is much smaller than that in flap. The ratio 
of the amplitudes is 

where u is defined in (B.20) and usually has a value of 0.1 to 0.2. 

(3) In hovering, if 0>8 C and k 3 r < ® / ^e l-*- m ^ c Y c l e ampli- 
tude is u 

A l - 1 



for lag motion and 



K3R 



for flap motion. 

The following are the cases for 0 near 0 c in forward flight. 
The value of f2_ has strong effect on the behavior of blade due to 
periodic forcing function and the periodic parametric excitation. 
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(4) If ^p— 1/2 , the amplitude of the flutter mode is governed 

by (31). By taking U 2 =l, the amplitude of the lag motion is 
1/2 

y A(t 2 ). In other words, when flutter occurs, the lag motion 
1/2 

grows as y The time scale x 2 for A is yip, i.e., A varies 

more rapidly in the time scale \p (~l/y) for higher y. The 
stability near the origin [A=0] is governed by the sign of 

K * & ( Kir (9*80 * \ K$| X - [ (74a) 

or in general, the stability of the singular points of (31) 
is governed by the sign of 

K = fc. (?„-».! Ifrf- ) n«b) 

where 

— Ki K-3A 0 ft® "+ t 2-} 

S s - Ms + * kX, <75a) 

in which A=A is a singular point of (31) . (Note A =0 is a 
o o 

singular point) . If k of (74a) is less than zero, the origin is 
stable and if is also less than zero, since (31) is homo- 
geneous, the flutter mode is unlikely to be excited. 

(5) If flp=l, the flutter amplitude A is governed by (37) or 
(39). For small value of y, (small as compared to 0— © c ) , i.e. y^=e 
or y=e 3/2 with 0 9 being of order one, (if y<0(e 3/2 ) for 6-9 =0(e) 
forward flight will have little effect on the flutter of the blade) , 
the lag amplitude is 

Mu) = A (/»*<!') 
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i.e., the lag motion grows as y^ 2 and the length of time for A 
to grow is proportional to 1/y 2 ^ 2 . For larger values of y (large 
as compared to 0-0 c ) , i.e., y^=l or with 02 being of order 

one or smaller, the lag amplitude becomes 

e ACM * r ft 


i.e., the lag motion grows as y and A varies more rapidly in the 
time ip (~l/y 2 ) for higher y. The origin (A=0) on the phase plane 
is no longer a singular point since (39) is inhomogeneous. There- 
fore, depending on the distance from the singular point to the 
origin, the stability prediction by the linearized theory such 
as (41) may not be meaningful. Let A=A o be a singular point of 
(39), then the stability of the singular depends on the sign of 
k [defined in (74b)] where 


Sp * ^ Kq+>*K„A. + KjA 0 


Since (39) is inhomogeneous, A q ^ 0 , i.e., the flutter mode is 
always excited. 

(6) If 2, the flutter amplitude is governed by (44). 
F 

For small values of y,i.e., y.^e^ 2 or y=e with 0 2 being of 
order one, the lag amplitude is 

e* Ate.) = ^ A (fit) 


1/2 

For large y, i.e., y or V =e 
the lag amplitude is 




wi th 


©2 being of order one, 
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The stability of the singular point depends on the sign of k 
[defined in (74b) ] in which 

K fe -+ i (ri f -2.) - iVu 

where A=A q is a singular point of (45) . Since A q =0 because (45) 
is homogeneous, the flutter mode is always excited. 

(7) If ft ±1/3 or 3, the flutter amplitude is governed by (50), 

F 

1/2 

or (54). Taking y^=l, i.e., y=e , the lag amplitude is 

£* KtO = / 


The stability equation for the singular point is also (74b) in 
which 




1 (75c) 

Ss- = i H, A 0 r ij'KyAa 


for ^ F “l/3=0(e) and A Q e 1VX 2 is a singular point of (50), 


K 2 ce-e c ) 4 t 


(75d) 


* K % A o 

for ft =3=0 (e) and A e 1VT 2 is a singular point of (54). In the 

F O 

latter case, (54) is inhomogeneous, the flutter mode will be 
excited. 


(8) Others. The amplitude equation (55) is similar to that 

of hovering. Therefore, the conclusion (3) also holds in this case 

2 

if k 2R^~®c^ re P^- ace ^ by K 2R (0-9 c )+y k 6r* 
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(9) When evaluating according to (B.65), if the quadratic 

terms in 3 12 and are ne 9 lected, and in (75) A q (or A Q ) are set 

to zero, k=0 of (74) gives the stability boundary of the linearized 
solution of (1) and (2) . 

(10) If < 3r <0, ftp is not near one, two, or three and if the 
blade is operated under the condition k<0, the flutter mode will 
not be excited. Therefore, in practical design, ftp=l, 2, or 3 
should be avoided. 

(11) The amplitude of the blade motion is of order e | A | , 
how e relates to the advance ratio y depends on the value of ftp. 

In other words, for the same value of y, the magnitude of the 
response can be of different order. In particular, if ftp=l, 2, or 
3, the flutter mode is always excited; one cannot just from the 
observation of the motion of the blade conclude whether the 
blade is in the region of stability or instability. Because the 
motion cannot simply be the forced response, for different value 
ft , the response is magnified to different magnitudes. 

Numerical results are presented in the figures. The 
induced flow is approximated by 



where a is the solidity ratio (Ref. 4) . Throughout the computation, 
the following parameters are used • 

- 0 . 0 ^ 

°-ol 

a - 
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For practical interest, 
2 

and 0^~1. 


we shall only compute the cases 1+0 2 ~1.1 

p 


In Figures 2 and 3, the locus of the critical conditions for 

various 6 c are presented; y=5 in Fig. 2 and y=10 in Fig. 3. These 

loci are similar to that of Fig. 4 of reference 4. The region 

within the locus is unstable according to linear aeroelasticity 

theory. For y=10, such regions are much larger than that for y = 5 

for the same 6 C . The numbers given on this locus are the values 

|k 2 r //|C 3rI ' for the corresponding parameters at 0=0 C . The solid 

line is a portion of the locus where < 3R <0 and the dotted line is 

That is, over the solid line portion of the locus, there 

exists a stable limit cycle of amplitude 2(6-0 O 1 / 2 Ik od /k, t J 1 ' /2 

C 2R 3R 1 

for lag motion and 2 (9”9 C ) ^ 2 | k 2r/ k 3rI l u l ^ or when 0>0 C * 

I u | is usually 0.1~0.2 because of low aerodynamic damping in lag. 

The blade is unconditionally stable if 6<0 . Over the dotted por- 

c 

tion of the locus, the blade is ' unconditionally unstable if 9>0 C and 
conditionally stable if 0<0 c , i.e., if the disturbance in |A| is 
smaller than I K 2R^ K 3R I ^ die out and if it is larger 

than I k 2r/ k 3rI ’ grow. In these two cases, the limit 

cycle amplitude can be quite large for 6>0 c , because |< 2 R /K 3R |a 
is quite large, i.e., the nonlinear effect is weak. For such a 
blade, it is not desirable to operate above the critical condition. 


Figures 4 and 5 are the lag response of the blade in hovering 
for a given initial disturbance, obtained by 

method of numerical integration. In both cases, j 1+^ 2 =1. 09544 , 
y=5, 0=0.3 and 6 c =0.25. 0^=1.0017 for Fig. 4 and 0=1. 19747 for 
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Fig. 5, while the numerical values for k^ r are -0.0074 and 0.001 
respectively for both cases. It is clear that the motion in Fig. 4 
approaches a limit cycle oscillation while the amplitude in Fig. 5 
grows steadily. 

The case considered in Fig. 4 has a flutter frequency ^ F =l. 00418 
It is expected that A is to be governed by (37). In this case, 

Hx - O.o(m+ 0.0l7ofcc 

-o.oo"f < H4' -t o.oioi t 

r-6.ou^\-0.00«*l i 

K<, _ -a o H W - »• oo 1 i 

* 

K, t0 = 1 

0.0* t»*l + O.OOlV* ' 

0. 0 I I b~! * <• 07-^1* < 

fl F _| = 0.OO<H8 

The time variation of |a| are plotted in Fig. 6. As predicted 
in Sec. IV. 2, when ft_~l, the amplitude of the response is of 

r 

for (0-0 ) /y^^ of order 1 and the time scale is proportional 
c 

to l/y^/^. Figures 7 and 8 are the lag responses of the blade 
obtained by numerical integration having the same physical parame- 
ters as that of Fig. 6. y=0.03 for Fig. 7 and y=0.04 for Fig. 8. 

(0-0 c ) /y^' 5 are held to be one. The lag motion is normalized 

1/3 2/3 

by dividing by y / , and \p is normalized by multiplying by y 7 . 

It can be seen that the amplitudes of the two curves are very 

2/3 

similar in the time scale x 2 =y rp. 


-36- 


REFERENCES 

1. Young, M.J., "A Theory of Rotor Blade Motion Stability In 
Flight," J.A.H.S., Vol. 9, No. 3, July 1964, pp . 12-15. 

2. Hohenemser, K.H. and Heaton, P.W., Jr., "Aeroelastic Insta- 
bility of Torsionally Rigid Helicopter Blades," J.A.H.S., 

Vol. 12, No. 2, April 1967, pp. 1-13. 

3. Tong, P., "Nonlinear Instability of Rotor Blades in Flap-Lag 

Using Hohenemser* s Equations," Unpublished Report, MIT Aero- 
elastic and Structures Research Laboratory, May 1971. 

4. Ormiston, R. A. and Hodges, D.H. , "Linear Flap-Lag Dynamics of 
Hingeless Helicopter Rotor Blades in Hover," J.A.H.S. 

Vol. 17, No. 2, April 1972, pp. 2-14. 

5. Cole, J., "Perturbation Method in Applied Mathematics," 
Blaisdell Publishing Company, Waltham, Massachusetts. 

6. Kaplun, S., "Fluid Mechanics and Singular Perturbations," 
edited by P.A. Lagerstrom, L.N . Howard, and C.-S. Liu, 

New York: Academic Press, 1967, 








40 




E' 


e- 


g M 

S § 


^ 3 
S E 


o 

o 

n 

o 

T3 

O 

A t— 


in 

rH 

o 

O 

• 

O 

o 

• 

n 

rH 

D 

n 


VJ' 

A 

c; 

in 


li 

«« 




-=r 

A 

in 

OS 

ON 

CM 

o 

• 

• 

o 

rH 

li 

U 

<D 

cm ca 


a 

VC 

+ 

o 



N 

55 

m 

o 

O 

M 

a 

H 

U' 

M 

ca 

Q 

a 


ii 

O 

kj* 

O 

cr 


n 


ca 

< 

c* 

o 


M 

+ 


p= 

M 

CM 

PC 

II 

O 

cd 


S uopow d^-[^ 


O'- 

5 uot^ow 


9 * 



41 - 


LJ 

”cr> 



w •» 


> H 


o o 


CQ 


< o 


It 


W 

O 

J T3 • 

o o 

t"- 


-=r 

< •» 

f- 

LTv 

On 

US O 

i— 1 

o 

• 

£h O 

rH 

H II 

II 

pL< b 

kJ* 


cz 

tt • 


O LTV 

* 

fe 11 


>“ 

-=r 

C5 

in 

!5 

ON 

m in 

o 

« nj 

• 

w 

rH 

> O 

H 

O ll 


W o 

icvj ca 


05 

-3- a od 

+ 

W O 

IrH 

fe 


25 


O 2 

m 

HO 

O 

2 EH H 

n 

KO E-» 

k> 

*rj s h 

ca 

H Q 

c; 

. U 2 

n 

rH < O 

k> 

2 ^ o 

cr 

c 

ll 

o Q J 

ca 

* H Z < 

er 

m 3 o 


C t— 1 

« 

g (X, ^ 

t= 

e < m 

C\J 

zi ^ « 

ll 

c fe ° 

cd 

o 


2 in 


• 


o 


M 





9 UOT^OW 


5 uoTq.ow 



- 42 - 





FIG. 7 NORMALIZED FLAP AND LAG MOTIONS VS. SLOW TIME SCALE 
T,(=y 2 / 3 iJ0 BY DIRECT NUMERICAL INTEGRATION OF 



FIG. 8 NORMALIZED FLAP AND LAG MOTIONS VS. SLOW TIME SCALE 

2/3 

t 2 (=vi ' <J>) BY DIRECT NUMERICAL INTEGRATION OF 
EQS. (1) AND (2) FOR y=0.04. 
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Appendix A. Derivation of Flap— Lag Equations 

In deriving the equations of motion of the blade in flap-lag, 

• 

g , £,and 0 are used for flap, lag s and pitch angle respectively. ( ) 
denotes differentiation with respect to ^ and ij;=ftt with ft being the 
rotating speed of the shaft itself. Two coordinate systems are 
used, namely, (x,y,z) rotating with the blade about the shaft 
and (X,Y,Z) rigidly attached to the blade. The base vectors 
of the two coordinate systems are related by (Fig. 1) : 


l = GO X — 5v/x.p> (0^9 T, — /W £ 


(A-l) 


0 =: OP 8 ^ ~ K 

c* 

k — ^ p) S ** ^ ■+ 

Let the distance of a point on the blade from its support be x, 
or R£ where R is the length of the blade. The position vector, 
the velocity, and acceleration of the point can be written, as 


I R i. CAj ^ ^ (i ^ 


M " 


SI [- jg.] 


(A- 3) 


a “ si 

*+£ k ^ (»+£ ) ^ 3 (A- 4 ) 

"t t (i — *5 a ^Q> ] | 

Let m be the mass per unit length, then the inertia moment about’ 
the z-axis (the hub) is 
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=■ XJ7. 2, ^ -2 0-+^) onp k") (A_5) 

4 X5Y T. |i> “t- 0-+ aX- ^>3 


where 

1= ( « s f*td^ 

J O 


(A- 6) 


The flow velocity relative to the blade is 

v ft = ^ e ft [ on 0+^)1 - -s^ (4-+^ >j.i 


(A- 7) 


where V IND is the induced velocity. To evaluate lift and drag, 
it is more convenient to express V A in terms of (X,Y,Z)- 
coordinates, which move together with the blades, i.e. 

Ya= v*x -V y t +v* k (a - 8) 

It is clear that is the velocity component in the axial direction 
along the blade, while V and V are in the cross section of the blade. 

y z 

Y ua = -f^j$ 

using (A-l) we find 

V Y ~ -X £ (^+'c ) $<^4> ^'9 4 a f v|4Y ) u > ej 

+ 1 + c ^>7^ > 3j’ (a- io) 

V 2 = - Si r{X 4 ,u [ (i|^ ) wf> «p 0 - Av- Ohs) s^9_j 

+ 1 ( ,+ r ) i*.v-9 + ^ oo 9 3 ^ 
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We shall assume that quasi-steady two-dimensional aerodynamic 
theory holds for each cross section, i.e., the aerodynamic load 
vector can be written as 


where p A is the air density, "a" lift slope, c chord length, 


profile drag coefficient 


V T = V Y 1 •+ V a K (A ' 12) 

is the effective angle of attack, i.e., (Pig. l.C): 

« - W' ^.= Ik (A-13) 

v Y v Y 

and N is a unit vector normal to y T and in the plane of the blade 
cross section defined by 


u = 

r-> 


i<)fr Vy 

1 $ I 1V Y | 


(A- 15) 


If V >0, i.e., the flow is from the trailing edge to the leading 

y 

edge, this is called reversed flow. In our present 

2 C d 

consideration, we assume that V <0, and a ,,<<1, __o<<l. A 

y err a 

substitution of (A-12) through (A-15) into (A-ll) yields 


L a = (L„V i -v(L^i-^) A js 


where 
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IU a =4 f A «c {[ + % M Vr] 

*l%lv,\V e -V x Vy>>s}a^ 
mv (i v 2 + % ' Vy| vy]co& 
-[%lV Y lV*-V*V<f]A^9i 

(Lj) A -iL*) A ^| 

The aerodynamic moment at the hub (z-axis) is 

m a = S’i'L A e^ 

The equations of motion of the blade are 

OS + (S) 'X +I “i ^ ^ = 0 

in flap and 


(A- 17) 


(A- 18) 


(A- 19) 


(IS* t! A )‘ iS *1 »5 (A-20) 


in lag, where 8 is the preconing angle, and w fi , w are the 
bending frequencies in the flap and in the lag direction, is 

due to the elastic coupling. (Note that, depending on the con- 
struction of the blade, w 0 , , io a „ can be a function of the 

P C PC 

pitch angle) . In (A- 19) and (A-20) , it can be shown that by 

2 C d 

neglecting small terms, such as o compared to one, 
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X TO 1 T Bs^+fAp*/- 1 
+ 8 1 

- (A <*>£ 0+5 y ^ <0r*\ ) <-] C .* e) 

■+ f „ (J Op A-O^ ) A~03«^ (A _ 21) 

^H I t^ A ).JJ= I rf'[s^- J C' + 5 )(> <Ht 3 

* |bc^9+[A(1+/^(^Oa^cJ6 + ^^3 

- ^A co@> U + S> ■* /* C 

- jjUC>{3> (l+^^^/UlAUpW^OD^^) 0+i )-4- 5t u 4 (<|»*£) 

^ 60 (<!+!:) 

— [u^ ^ 0 + ^ ) 3 *- fy+'S ) * * 

where 

^ _ Jkac R 4 

• X 

is the Lock number 
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A*+ !.**<** 

e>=4 ( A 1 £ <A\ (A-22) 

Jo 

C = 4^d2, 

To further simplify (A-19) and (A-20 ) , we shall assume that 

elastic coupling w <0(e), the pitch angle 6, the static equilibrium 

p C“ 

flap angle B Q , and the static equilibrium lag angle 5 q are small, 

2 2 

so that 0 , 0 0 and Cq can be neglected as compared to 1. The 
quantities 3 Q and t 0 satisfy 

il >. + ($>.-£,0 + ^ -*) s ? 

r (A-23) 

We shall define the following non-dimensional parameters 

Si\ - ***/& 

d* = F 

&^Ab-2^>) <a-24> 

- 2 . = ^ (20 -A} 

T = f (6-2A) 



Substituting (A-23) into (A-19) and (A-20) , replacing 6 by 
B+Bg and z, by C+? 0 f i.e. using B» C to denote the perturbed 
motion flap and lag, respectively, and retaining only proper 
linear terms we obtain Eqs. (1) and (2). 
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Appendix B. Definition of Various Functions and Quantities 

Uf,^ = + t 'B- 11 

M(^5) = l-^a 5 t+-ns^+%p»- Y(i <b-2> 

■^,0^6)= ^. 0 (b-3) 

•S’cv.e) ~ lB . 3) , 

V^®) ^*ce)-sUi|/4 (b— 4 ) 

<V,0V,e} = (b-4) ' 

■$ J ('l',e) = g -©«»*</’) (B-5) 

■^OKe) = 4»cv,6'i /'f (b-5) * 

i 2 ct&) = +((£>-; (B_6) 

-ft’ u,e) = ft 1 (‘V,e)/‘} tB-6) • 

+|p.p-C9S] w5 +} 

-f 4 c 1 ^ (B-9) 

■fc+OM) =2.P>f> + e) •+|^ 1 


(B-10 ) 



(B-12 ) 





(leis-lpi-fis) o©<(- J 

l s (to; ^>gU ( f 1 *!p.B + Ip.»?5]5-')' 

£ ?*&>£- [*(£•«• + e ^,], w <|' 1 



W*a) - 


1 


/V/ 


c*]L -Rp^ A 


l+Rp- 


3>(*,*) 


(B-13) 


(B-14 ) 
(B-15 ) 


(B-16 ) 


(B-17) 


(B-18 ) 
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"DW — ^ 4 ^3x^0" °Y^) (b-19) 

4> SLu 


a = 




(B-20) 


IX* - — Cdwplex Cd>rtj*<joite o ^ K. 

Kp" Sip (A ^ (^6~ )"^p ^ ^ 

2 C Sip a 1 - (*$,- * - jSl f * 


T. \ * 


°V~ ( r^% ^ U ^ u ) + -J -2^p u*"lt 

H --(^JVVV-^o 
R = - (**% 4 3*0 ^Yo ^ 

P3= ~ A $(b2.Sif a + t5ipi 2 -J7psa. 

H- - C -<•(<- S?f Y>.-il^z) a. 


(B-21) 


(B-22) 


(B-23) 


^ - -2 t,#p(u*- g^Slp- yj£p Cu*4u) (B-24) 


(B-25) 


(B-26 ) 


(B-27) 


(B-28) 


(B-29 ) 
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V(,— 2. i- u <.-R.p^> ax + 2 J P<r^®) u ■% V»| 

4 (4^o- j9) Z Sp 4 £ j£ U^ M 4 2 U X J2p 

?s * a [cl C t **) * ■ + 1 Cc -**} 



(B-30) 


(B-31) 


(B-32) 


(B-33) 


(B-34 ) 


(B-35 ) 


(B-36 ) 
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(I 

+ [--^i 4 2 - f- ? F ^'*> „ 

P w =f [Tx'V| eu %Hc-fe)un| W -i%^|^ P i-%ej 
-» [a «% i - 2 u*§„ 2 F £ ^ &- f eX^Ta* «f 


(B 


4 rr j ( 

4 




(B 


p >^ f [-! 

+ i ia ^ l ***> 

'Rj — g [ f ®$|C "3 Pii" 1- (?®* j)?,,- jP'r j PoV i J-* 0 

~ a [-3 3®ft*. + ( c 'l®)p,t 4 3®o?,»i.- 5#^5i»i "2 Ce T,i.j«<, 

■ - I p 1 -' 

'la- ^ rti 


-|o)^, - r i - 


-37) 


-38) 


-39) 


-40) 


-41) 


(B-42) 
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TV = g [-§ 

(B-43) 

- f (Pu? u "^n X, y j. |se^ 

{B “ 44) 

- % ttM,4K + (& v^)*. 

-*• i Mi * 



( B— 45) 


(B-46 ) 


( B— 47) 


(B-48) 
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Tn = f if 9( V') 5* - 1 C»rOpf r |fc,< i +l|*-f )£- 

-fiiJV'K4 - fecvOW* I 

-t$®-t) V 1 f.lvOi +[^0 Vnt*(V &) 

4 (f 8- fu>% 

- I 1 ( &»' [C^T'O 1 f'l/^5 + jj* (B-49) 

-2. ( f.- |e)[CR T t|)3*^ -(Jl r ,)s u ^j* 0 - 2 cft f ^ 

4 £M<^ vrf w 

^ 0(3 Ph- | 8 r H )*^r!)^. (^I,„)+jp,M. 

-4 ( (t - f e>)l u (<; $„)*(£' (VI C;,t-r n ) _$*. ^ j 

* I b W% (>‘-V)-^-V] •* I u i [fa .- (£ -► M $ * )~\ (B_50) 

4 2t J2p (&x ^ M-') Jf <. + 2£%a^ 46 

■^[evo 

4 a^4,X)*; 
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- 1 U* ; 4 « r i (-e «*)$ ] . Ji 
■^-^F fi^n^+SZp «„ 

- (*6 - c) u + t - l&oU’^F u*-8f +f 0u*.;% 

^l“*e»«v>M£ 2M»] * ^ [-* T « + iw 

■* ^ ft»"$ &) [o - %) ( ‘^W&*V* * C + 3,1 )4?o] 57 F < B - 5 2 > 

-* | j> [ u *T^ 3 (>-%)r J ft, ^ 

-*• L 2 % ( «*3tt--»Fph)>. 

P», * II le vaipfe -aMf®-s)S* -f fas.i »* 

4 *t c -5 £ 3,, ft, s 

* |<&*V+ fcrV>| + M,t ^ 4 


(B-53) 



(B-54) 


= ? H - feM(c-ie^l) ft*- (Ki + f 

- Po(Mo 

+i eiv^ -» (p-*+ %)v*»} + [ 

t fc-'M (i?2>- S I W (6-f^j 

' +olt for fou*+ =Po^-^%) - n. f ^ 

F ) W" s -2? - f p^iO-^FJ 

"* f*» [' 4® -%) * C + 5 > C<- *{$l * 

^ ^ f . <** - ♦ «*% + a i $*««■(» . ; ?F ) 

- . ii[u(i-n F ^ t p 2> -xz f ^ j 

— ^.flf t^>„- |- ) j(W?p)<J z .j "(''^if)^ l . F j - ^ - ( u-> U*J (. 


(B-55 ) 


(B-56) 


(B-57) 
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Pib = ai[“f. 1 -i^(^1 + (i f,„- - (i-T-f )u^ 

+ (WJI F )(\^ - (l-J)f) (^ 3 .^-. £Slp 1 (B-58) 

"*■ yftc- 5;P)(u-tV + )t "t 


Tz r f [(&- »& J *t ) V -( i*~ i ' ; fc‘)s.r -f IV?<sJ 

■+ ( c - 1®) (V Cc * fa) (U <■ f (!>£« i- C 4$ir\ 


?, r f^CM) 

P J( , - \i; 

i. $p X ” -H^5 c 
^ l+-^p£. -flf +<ftf 3pC. 

H, » vrP, ♦?* 

Kj = - (JP^-P^Vk, 

K - )/k, 

>1,J= - ^ p* * ?£)/*, 


(B-61) 

(B-62) 

(B-63) 

(B-64) 

J=-0,4-,r (B-65) 

v»2- (B-66) 


(B-67) 



